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Abstract 

We construct a cosmological model with non-minimally coupled scalar field on 
the brane, where Gauss-Bonnet and Induced Gravity effects are taken into account. 
This model has 5D character at both high and low energy limits but reduces to 
AD gravity in intermediate scales. While induced gravity is a manifestation of 
the IR limit of the model, Gauss-Bonnet term and non-minimal coupling of scalar 
field and induced gravity are essentially related to UV limit of the scenario. We 
study cosmological implications of this scenario focusing on the late-time behavior 
of the solutions. In this setup, non-minimal coupling plays the role of an additional 
fine-tuning parameter that controls the initial density of predicted finite density 
big bang. Also, non-minimal coupling has important implication on the bouncing 
nature of the solutions. 
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1 Introduction 



The idea that we and all standard matter live on a brane embedded in a higher dimen- 
sional bulk has attracted a lot of attention [1-3]. In this view point, extra dimensions are 
accessible only for graviton and possibly non-standard matter. The setup of Randall and 
Sundrum (RSII) considers the observable universe as a 3-brane with positive tension em- 
bedded in five dimensional anti-de sitter bulk. In the low energy limit, the 5D graviton is 
localized on the brane due to the warped geometry of the bulk. The notion of AdS / CFT 
correspondence help us to understand this property since the 4D gravity is coupled to a 
conformal field in the RS model [4,5]. 

The effect of the bulk on the brane can be determined by the effective mass Ai of the bulk 
fluid that is measured by a bulk observer at the brane. For spherically symmetric brane, 
this mass can be considered as effective gravitational mass of the bulk. This mass depends 
on the brane scale factor and the proper time on the brane. In the case that the bulk 
observer is comoving with the bulk fluid, the mass is assumed to be comoving. However, 
for matter components such as a bulk radiation fluid, there is no comoving observer [6, 7]. 

On the other hand, the model proposed by Dvali, Gabadadze and Porrati (DGP) 
is a radiative correction, the bulk is a flat Minkowski spacetime, but a reduced gravity 
term appears on the brane without tension. This model is different in this respect that 
it predicts deviations from the standard 4-dimensional gravity over large distances. In 
this scenario, the transition between four and higher-dimensional gravitational potentials 
arises due to the presence of both the brane and bulk Einstein terms in the action. Exis- 
tence of a higher dimensional embedding space allows for the existence of bulk or brane 
matter which can certainly influence the cosmological evolution on the brane. This model 
has a rich phenomenology discussed in [8]. Maeda, Mizuno and Torii have constructed 
a braneworld scenario which combines the Randall-Sundrum II ( RSII) model and DGP 
model[9]. In this combination, an induced curvature term appears on the brane in the 
RSII model. This model has been called warped DGP braneworld in literature [10]. The 
existence of induced gravity term leads to a self-accelerating branch in the brane evolu- 
tional,^]. 

Braneworld model with scalar field minimally or non-minimally coupled to gravity have 
been studied extensively(see[13] and references therein). The introduction of non-minimal 
coupling (NMC) is not just a matter of taste; it is forced upon us in many situations of 
physical and cosmological interests. For instance, NMC arises at the quantum level when 
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quantum corrections to the scalar field theory are considered. Even if for the classical, 
unperturbed theory this NMC vanishes, it is necessary for the renormalizability of the 
scalar field theory in curved space. In most theories used to describe inflationary scenar- 
ios, it turns out that a non-vanishing value of the coupling constant cannot be avoided. 
In general relativity, and in all other metric theories of gravity in which the scalar field is 
not part of the gravitational sector, the coupling constant necessarily assumes the value 
of | . The study of the asymptotically free theories in an external gravitational field with 
a Gauss-Bonnet term shows a scale dependent coupling parameter. Asymptotically free 
grand unified theories have a non-minimal coupling depending on a renormalization group 
parameter that converges to the value of | or to any other initial conditions depending 
on the gauge group and on the matter content of the theory. An exact renormalization 
group study of the A0 4 theory shows that NMC= | is a stable infrared fixed point. Also 
in the large N limit of the Nambu-Jona-Lasinio model, we have NMC= |. In the O(N)- 
symmetric model with V = A0 4 , NMC is generally nonzero and depends on the coupling 
constants of the individual bosonic components. Higgs fields in the standard model have 
NMC= or |. Only a few investigations produce zero value(for a more complete discus- 
sion of these issues we refer to papers by V. Faraoni, specially Ref. [14] and references 
therein). In view of the above results, it is then natural to incorporate an explicit NMC 
between scalar field and induced Ricci scalar on the brane. 

On the other hand, in a braneworld scenario, the radiative corrections in the bulk 
lead to higher curvature terms. At high energies, the Einstein-Hilbert action will acquire 
quantum corrections. The Gauss-Bonnet (GB) combination arises as the leading bulk 
correction in the case of the heterotic string theory [15]. This term leads to second-order 
gravitational field equations linear in the second derivatives in the bulk metric which is 
ghost free [16-18], the property of curvature invariant of the Gauss-Bonnet term. 
Inclusion of Gauss-Bonnet term in the action results in a variety of novel phenomena 
which certainly affects the cosmological consequences of these generalized braneworld 
setup, although these corrections are smaller than the usual Einstein-Hilbert terms [19- 
21]. Moreover, the zero mode of graviton has been localized in the GB model [22]. The 
cosmological evolution corresponding to RS model in the presence of a bulk GB term has 
been considered in [17,23-28] see also [29]. Also the case of minimally coupled scalar field 
with GB gravity has been discussed extensively [30-33]. 

In the presence of GB term with induced gravity, there are different cosmological scenarios, 
even if there isn't any matter in the bulk[34]. In this paper, we generalize the previous 
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studies to the case that a scalar field non-minimally coupled to induced curvature is present 
on the brane in the presence of radiative corrections. We first review briefly the setup of 
Brown-Maartens-Papantonopoulos-Zamarias(BMPZ)[35]. Then we generalize this setup 
to the more general framework of scalar-tensor theories. We show that relative to BMPZ 
scenario, there are several interesting features which affect certainly the cosmological 
dynamics on the brane. Since Gauss-Bonnet and induced gravity effects are related to 
two extremes of the scenario (UV and IR limits), inclusion of sringy effects via Gauss- 
Bonnet term leads to a finite density big bang [35]. This interesting feature has been 
explained in a fascinating manner by T-duality of string theory [36]. On the other hand, 
non- minimal coupling itself accounts for a non-singular soft big bang scenario [37]. In 
our setup, existence of non-minimal coupling of scalar field and induced gravity on the 
brane, controls the initial density of this finite big bang scenario. In other words, in this 
framework, non-minimal coupling with its special fine-tuning ( see [38] and references 
therein), plays the role of a parameter that can control density of matter fields at the 
beginning of the universe. As we will show, incorporation of both GB and non-minimal 
coupling effects will enhance special characters of BMPZ scenario. 



2 DGP Inspired Scalar- Tensor Theories 



The action of the DGP scenario in the presence of a non-minimally coupled scalar field 
on the brane can be written as follows [39] 



2k 2 J 



d 5 x 



,(5) 



2A, 



+ 
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J d A x^{a((t))R - 2/t^V M 0V i ,0 - AkIV{4)) - Ak\\ 



(1) 



2/=0 



where we have included a general non-minimal coupling a(0) in the brane part of the 
action, y is coordinate of the fifth dimension and we assume the brane is located at y — 0. 
g^AB is fi ve dimensional bulk metric with Ricci scalar R^\ while is induced metric on 
the brane with induced Ricci scalar R. qab and g^ are related via g^ = 5^ A 5 V E ' qab- A 
is the brane tension (constant energy density) and r is the cross-over scale that is defined 
as follows 

..2 i rl 

(2) 
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The generalized cosmological dynamics of this setup is given by the following Friedmann 
equation [34,40] 



where e — ±1 is corresponding to two possible branches of DGP cosmology and T is the 
bulk black hole mass which is related to the bulk Weyl tensor. This mass, as generalized 
dark radiation, induces mirage effects in the evolution and the gravitational effect of the 
bulk matter on the brane evolution can be described in terms of this mass as measured 
by a bulk observer at the location of the brane (the DGP limit has a Minkowski bulk 
A 5 = with T = 0). A part of the effects of non-minimal coupling of scalar field with 
gravity is hidden in the definition of the effective energy density. Assuming the following 
line element 

ds 2 = q^dx^dx" + b 2 (y, t)dy 2 = —n 2 (y, t)dt 2 + a 2 (y, t)^ijdx l dx^ + b 2 (y, t)dy 2 , 

where 7^ is a maximally symmetric 3-dimensional metric defined as 7^ = 5ij + k , 
the energy density of non-minimally coupled scalar field on the brane is given as follows 

[39,41] 

P<t> 



-<P 2 + n 2 V{<f)) - 6a'H<p 
2 



(4) 

y=0 



where H = ^ is Hubble parameter, a' — % and = 

If we consider a flat brane (K=0) with A = and also a Minkowski bulk (A 5 = 0, 
T = 0), then we can write equation (3) as follows 

where p 0( p = \4> 2 + n 2 V(4>) . The DGP model has two branches, i.e e = ±1 cor- 
responding to two different embedding of the brane in the bulk. The behavior of two 
branches at high energies and low energies are summarized as follows: 
In the high energy limit we find 

DGP{±) - H2 = ^» + ^ (6) 

while in low energy limit we have 



DGP{+) : H — ► — — r - 2 



ra{4>) ct(4>) 
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DGP(-): H = 0. 
In terms of dimensionless variables introduced in [35] 



rni , rni rni t 

h = Hr, p = —p, p = —po<p, o- = —A, r = -, 
6 6 6 r 



we find 



(7) 

(8) 
(9) 

~2a(4>) a(4>)dT ' 2a(0) ' ^ 

Here the negative root is not suitable since in the limit of p + p! — ► 0, with this sign 
one cannot recover the low energy limit of the model highlighted in (7). 



h 2 



h (p + p') 2hn\ da 



a((f)) cx((f>) ot(4>) dr 
The solutions of this equation for h are as follows 



h = ±- 



k\ da yfl T Ak\% + A{k\%)* + 4a(0)(/i + p') 



+ 
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Figure 1: Two possible branches of DGP-inspired non- minimal model. The non-minimal 
coupling of scalar field is assumed to be positive and the brane is considered to be tensionless, 

(7 = 0. 

Figure 1 shows the behavior of these solutions with some specific values of non-minimal 
coupling 3 . The upper sign in relation (10) is related to DGP(+) and lower sign for DGP(- 
). It is seen that there is a late-time self-acceleration in the DGP(+) branch similar to 

3 To plot all of the figures in this paper, equation (41) acts as a condition on the values that a 
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Figure 2: Two possible branches of DGP-inspired non- minimal model. The non- minimal 
coupling of scalar field is assumed to be negative and the brane is considered to be tensionless, 

(7 = 0. 



the minimal case, however, in minimal case when /i — > then h — > 1, whereas in the 
non-minimal case, in this limit i.e + — > 0, we have h — > 7^pj~ I 11 DGP(+) 
branch, the endpoint is a vacuum de Sitter state and Anti de Sitter state for positive and 
negative non-minimal coupling respectively whereas in the DGP(-) branch, the endpoint 
is a Minkowski state. The effect of non-minimal coupling in this case is to shift the end 
point of DGP(+) branch. Depending on the value that a(<p) can attain [38], the late-time 
acceleration of the universe can be fine-tuned properly. 

Now for A 7^ 0, the solutions of dimensionless Friedmann equation are as follows 

, = 1 4 da y/l?4gg + 4(«ffi)2 + 4a(flfc + ;*' + <7) 

2a(</>) a{4>) dr 2a(0) ' 1 j 

These solutions are shown in figures 3 and 4 with c^O. For negative tension in DGP(+) 
branch, the endpoint is a vacuum de Sitter state and there is a self-acceleration whereas 
in DGP(-) branch the solutions terminate at finite density (fig. 3). For positive tension, 
both of the solutions (DGP(±)) have self-acceleration and the endpoints are vacuum de 
Sitter state (fig.4). 

The existence of the energy density (A) on the brane gives rise to a shift of the solutions. 

can attains. The possible values of 7 extracted from observational data are shown in table (1). Using 
SNIa+LSS-|-H(z) test, we obtain a > 0.01 and a < —0.01 approximately. For simplicity in drawing 
figures we have assumed that scalar field has no dynamics, i.e. ^ = 0. The time dependent non-minimal 
coupling will be discussed at the end of the paper. 
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Moreover, existence of a(<f>) leads to further shift of these solutions. 





10 20 30 40 50 60 70 80 90 100 
energy density (n+n') 



Figure 3: Two possible branches of DGP inspired non-minimal model with negative tension 
brane and positive non-minimal coupling. We have set a(4>) = 0.01,0.02 and a = —4. 



DGP model is an IR modification of general relativity. In the UV limit, stringy effects will 
play important role. In this viewpoint, to discuss both UV and IR limit of the scenario 
simultaneously, the DGP model is not sufficient and we should incorporate stringy effects 
via inclusion of the Gauss-Bonnet terms. 

3 Gauss-Bonnet Braneworlds 

The Gauss-Bonnet term with coupling constant (3 is written as follows 

L GB = - AR^R {5)ab + R { l d R {5)abcd 

where is the curvature scalar of the 5-dimensional bulk spacetime. These corrections 
have origin on stringy effects and the most general action should involve both Gauss- 
Bonnet and the Einstein-Hilbert term in 5D theory. The GB term is present only in the 
bulk action 



S bulk = ± J d'xf^) \r^ - 2A 5 + (3(R^ 2 - 4R%R^ b + R^R^) 



(12) 
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Figure 4: Two possible branches of DGP inspired non-minimal model with positive brane 
tension and positive non-minimal coupling. We have set a(<j>) = 0.01,0.02 and a = 10. 

where (5 is the Gauss-Bonnet coupling which can be positive or negative in the classical 
GB theory. If (3 is negative, it has been seen in [42] that this braneworld model leads 
to antigravity or tachyon modes on the brane. However, in the presence of a bulk scalar 
field, these effects are not present even with negative (3. 

The Friedmann equation in the presence of Gauss-Bonnet effects is as follows [43] 



H 1 



C+ + C_ - 2 



8[3 



K 



where 



C ± = 



(13) 



(14) 



2 'V 2 

This equation is a cubic equation with three possible roots. For p > there is only one 
real root. 

The behavior of GB model at high and low energies are as follows [35] 



and 



H > a 5 _> H 2 oc p3, 
H < oT^ -> H 2 oc p 2 



high energy limit 
low energy limit. 



(15) 
(16) 
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4 Gauss-Bonnet Induced Gravity with Non- Minimally 
Coupled Scalar Field on the Brane 



As we have explained, Gauss-Bonnet effect is a high energy stringy effect. On the other 
hand, non-minimal coupling of scalar field and induced gravity on the brane is forced 
upon us from several compelling reasons. Some of these reasons have their origin on 
pure quantum field theoretical considerations [14]. Then it is natural to incorporate both 
Gauss-Bonnet and non-minimal coupling effects to have a more reliable framework for 
treating cosmological dynamics. The action of the GBIG (Gauss-Bonnet term in the bulk 
and the Induced Gravity term on the brane) scenario in the presence of a non-minimally 
coupled scalar field on the brane can be written as follows 



S = — d 5 x 



24 



,(5) 



- 2A 5 + p(R^ 2 - AR { ^R^ )ah + R^R®^ 



+ 



(17) 



y=0 



where (3 and r are the GB coupling constant and IG cross-over scale respectively. The 
relation for energy conservation on the brane is as follows 



p + p 4> + 3H(1 + oo){p + P4> ) = 6a'<j)(H 2 + ^) 



(18) 



where uj = with p and p pressure and density of ordinary matter. Since for ordinary 
matter, p + 3H(p + P) = 0, the non-minimal coupling of the scalar field and induced 
curvature on the brane leads to the non-conservation of the scalar field effective energy 
density [41]. 

The cosmological dynamics of the model is given by the following generalized Friedmann 
equation 



K k 2 

ra(0)/J 2 + m(0)— --|(p + p + A) 



1 2 



• (19) 



This equation describes the cosmological evolution on the brane with tension and a non- 
minimally coupled scalar field on the brane. The bulk contains a black hole mass and a 
cosmological constant. \1/ is defined as follows 

A 5 T 



^ + 2(3^ z = — - + — . 



(20) 
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If (3 — 0, the model reduces to DGP model, while for r = we recover the Gauss-Bonnet 
model. Here we restrict our study to the case where bulk black hole mass vanishes, 
(T = 0) and therefore \I> + 2[3fy 2 = The bulk cosmological constant in the presence 
of GB term is given by A 5 = —jj + -jf, where I is the bulk curvature. For a spatially flat 
brane (K — 0), the Friedmann equation is given by 



ra{<P)H 2 --A( p + P0 + A) 



3' V 2 

We define the following dimensionless quantities 

8(3 _ r 2 

where the dimensionless Friedmann equation takes the following form 

da{4>) 



1 + 7 (V + ^) {h 2 -i) = a(<P)h 2 -(fi + fi l + a-2 



(21) 



(22) 



(23) 



To find cosmological dynamics of our model, we should solve this equation in an ap- 
propriate parameter space. In which follows, we consider Minkowski and AdS bulk and 
investigate their cosmological consequences. 



4.1 Minkowski Bulk(^ = 0) with Tensionless Brane(cr = 0) 

In this case, the non-minimal GBIG Friedmann equation takes the following form 



(l + 7 fc 2 ) V = \a(<P)h 2 -(// + //- 2^-hn 2 )] 2 . 



It is straightforward to show that in this case 
d(n + n') (l + 7 /i 2 )(3 7 /i 2 + l) 



d(h 2 ) 2 («(0)/ i 2 -( /U + /i') + 2^K) 



a(<P) 2 h 2 + M<P) dj ^hK 2 - + //)(«(«/>) + 1^1) + 2{^n 2 f 



2 (c<0)/i 2 -(/x + //)+2^^ 2 ) 

The initial Hubble rate and density are given by substituting the result of - 
to equation (24). This leads us to the following relations 



(24) 



hi 



a{<f>) + Va(0) 2 - 3 7 + 6 7 ^«i 
3 7 



(25) 



in 



(26) 
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2a(<Pf - 9a(0) 7 + 18a(0)^/#y + 2(a(0) 2 - 3 7 + 67 



da(cj>) „2\ 2 



27 7 2 



(27) 



Before proceeding further, we should stress on two important points here: firstly, the 
presence of GB term removes the big bang singularity in this setup, and the universe starts 
with an initial finite density. Gauss-Bonnet effect is essentially a string-inspired effect in 
the bulk which its combination with pure DGP scenario leads to a finite big bang proposal 
on the brane. A consequence of string inspired field theories is the existence of minimal 
observable length of the order of Planck length [44-46]. One cannot probe distances smaller 
than this fundamental length. In fact a string cannot live on the scale smaller than its 
length. This feature leads us to generalize the standard Heisnberg uncertainty relation to 
incorporate this Planck scale effect [47,48]. The existence of this minimal observable length 
essentially removes spacetime singularity and acts as a UV cutoff of the corresponding 
field theory( see for instance [49] which discusses inflation with minimum length cutoff. 
See also [50]). So, in principle existence of a finite density big bang is supported at least 
in this viewpoint [51], see also [52]. Secondly, non- minimal coupling of scalar field with 
induced gravity on the brane controls the value of the initial density. This is not the 
only importance of non-minimal coupling of scalar field and induced gravity. In fact non- 
minimal coupling provides a mechanism for generating spontaneous symmetry breaking 
at Planck scale on the brane [53]. In this respect and based on the arguments presented at 
the introduction on the importance of the non-minimal coupling, non-minimal coupling 
of scalar field and induced gravity on the brane itself is a high energy correction of the 
theory and it is natural to expect that this effect couples with stringy effects in Planck 
scale. In fact in this setup, we encounter a smoother behavior due to Gauss-Bonnet term 
( a finite density big bang) and the late-time effects of non-minimally coupled scalar field 
component. These effects together provide a more reliable cosmological scenario. 

The behavior of h with respect to //+// is shown in figure 5. This figure shows also that 
the GBIG1 and GBIG2 branches have self-acceleration for some positive values of non- 
minimal coupling in the same manner as the DGP(+) branch, whereas GBIG3 branch 
similar to DGP(-) has no self-acceleration. This is similar to pure DGP or GB model 
alone where there is a big bang singularity. The self-accelerating GBIG2 branch is not 
a physical solution since it is accelerating throughout its evolution. For negative values 
of non-minimal coupling, the GBIG3 and GBIG2 have self-acceleration while GBIG1 has 
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not such a property. Now it is easy to show that 

(/x + //)—> oo : 7^0 



: 



7 



4(1 -2^4) 



In the minimal case, the maximum value of 7 leads to a minimum value of Here, 
in the presence of non-minimal coupling of scalar field and induced gravity with 7 ma x = 
4fi-2dL 2) ) we cannot conclude that this leads to a minimum value for h. Using equation 
(26), the value of hi for 7 max is given as follows 



2(l-2f«2) 



«(0) 



(28) 



When 7 = 7 max and (// + //) = 0, there is a vacuum brane with de Sitter expansion. The 
h asymptotic value is obtained when /i + y! — > in equation (24) 



"'00 ~ 



(27 ~« 2 ) h 4 



-ht + 



hi 



0. 



(29) 



This equation has four non-zero roots that two of them are negative and unacceptable. 
When 7 — ► 0, we should have the non-minimal DGP model. 












GBIG1 (a=0.2) 






'•^-GBIG2 (a=-0.2) 














GBIG3 (a.0.2) 




: GBIG1 (cfc-0.2) 










12 3 4 5 
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Figure 5: Solutions of the Friedmann equation for a tensionless brane in a Minkowski bulk. For 
clarification, we have plotted the figure in two different scales to highlight intermediate points. 
Also we have set 7 = 0.003 and a{4>) = 0.2, -0.2. 



From equation (24) one can deduce 



da 



(/i + n') = a((j))h 2 - h^h 2 + 1) + 2—Klh, 

dr 



(30) 
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where < h < hi. Since (/i + /i')oo=0, from equation (30) it follows that 

7 = — • (31) 



hi 



By expanding // + // to first order in /r — ft, 2 ^, we find 

2(^(0)^ + 2f^ 



/> 2 = ^ + 



a(0)/ ioo (2 - 6f «| - a(0)/ioo) - 8(f ^) 2 + 4g«| 



(// + //). (32) 



In comparison with equation (5), we find the following effective 4-dimensional Newton's 
constant 



G = 



a(4>)h 



, nda 2u 



a(0)/i oo (2 - 6f «2 - a(0)/ioo) - 8(f k 2 ) 2 + 4f « 2 



a(0)G 5 



(33) 



where G 5 = ^ and C = ^ are five and four dimensional gravitational constant respec- 
tively. From equation (33) we obtain a relation between Mf and M 2 as follows 



(a(0)r 2 tf 2 + 2f k\tH q ) 



_a(<P)rH (2 - 6f k 2 - a(<j>)rH ) - 8(f^ 2 ) 



2 I a da 2 



a(0)M 2 



(34) 



Here ifoo ~ ifo, and we see the important role played by non-minimal coupling in this 
setup. In principle, one can fine tune the value of the non-minimal coupling such that 
fundamental scale of the bulk be reduced to values in the range accessible for next gen- 
eration of accelerators. To compare with DGP(+) limit, when ji + // — > 0, that is at 



late-time, we have rHo 
to Ml ~ ^ 

o r 



1 

a(4>) 



da 



' a((f>) dr 



, therefore this equation in this limit reduces 



This is an interesting result since there is no effect of non-minimal cou- 
pling in the non-minimal DGP(+) limit at late-time and the relation between M| and 
M 2 is the same as minimal case. This is not surprising since essentially a part of motiva- 
tion for inclusion of non-minimal coupling has its origin on the quantum field theoretical 
considerations( the renormalizability of quantum field theory in curved background and 
quantum corrections to the scalar field theory). In this view point, non-minimal coupling 
shows its importance mainly in the high energy UV sector of the theory while apparently 
DGP(+) gives IR sector of the theory free of stringy and strong quantum field theoretical 
effects. We should stress that non-minimal coupling of scalar field and induced gravity on 
DGP brane modifies cross over scale[39]. The above argument is restricted to the limit 
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// + //—> which is related to the late time stage of evolution. In equation (34), when 
«((/>) attains different values, M 5 increases or decreases relative to its value in DGP limit. 

2(1— 2 — k 2 ) 

When rH — > a ^ 4 , M 5 increases. We should stress that these results are sensitive 
to the sign of the non-minimal coupling explicitly. 



4.2 Minkowski Bulk (if) = 0) with Brane Tension(cr ^ 0) 

In this case the Friedmann equation is as follows 



[l + 1 h 2 ) 2 h 2 



a 



(0)/> 2 -(/i + /i' + a-2^/^) 



(35) 



The effect of brane tension is similar to considering a cosmological constant on the brane. 
In this case, there are three possible solutions (GBIG 1-3). To find initial and final Hubble 
rates and density, we set ^ = 0. The initial and final Hubble rates denoted by hi 
and h e respectively, are given by 



q(0) ± yW) 2 - 37 + 6 7 ^j 
37 



(36) 



and the initial and final density are calculated as follows 

_ 2a(0) 3 - 9a(0) 7 + 18a(0)^ 2 7 ± 2(a(0) 2 - 3 7 + 6 7 ^ 2 ) 



(fJL+fJL') 



27 7 2 



(37) 

where the plus sign is for initial state and minus sign shows the final state. These points 
have b! = 0. The point h = is the place which GBIG3 loiters, that it is given by 



(/i + fi')i = -a. 



(38) 



and is given by 



7 2 a oo T 



1 + 2aa - 4(f K 2 ) 2 



hi 



r 



(39) 



to obtain this relation we have set + = in equation (35). 

The relation between /i and // + // is shown in figure 6 for a negative tension brane. 
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In this figure for positive non-minimal coupling, the GBIG1 and GBIG2 branches start 
with an initial Hubble rate and density (hi, (// + whereas the GBIG3 branch doesn't 
remove the big bang singularity. GBIG1 and GBIG3 terminate at a finite Hubble rate 
and density (h e , (jj + /i')e), whereas GBIG2 terminates in a vacuum de Sitter state. On 
the other hand, for negative non-minimal coupling, the GBIG2 branch has a big bang 
singularity and this is a self-accelerating solution. GBIG1 and GBIG3 start without 
big bang singularity (hi, (fi + //)«) and both of them have self-acceleration but GBIG3 
throughout evolution loiters and then evolves. A universe which undergoes a period 
of loitering is an attractive alternative to standard cosmologies. Generally a loitering 
universe is an expanding Friedmann universe that undergoes a phase of slow expansion 
with redshift of (z ~ 3 — 5). It is believed that the large scale structure of the universe is 
formed during this semi-static phase [54]. 




10 20 30 40 50 60 70 80 90 100 5 10 15 

energy density (li+li r ) energy density (n+n') 



Figure 6: Solutions of the Friedmann equation with a negative tension brane in a Minkowski 
bulk. We have assumed 7 = 0.003, a = —4, ip = 0, a{4>) = 0.2 and —0.2 for solid and dot curves 
respectively. The figures are plotted with different scale to highlight important features. 
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Figure 7 is the result for a positive tension brane. Here with positive non-minimal 
coupling, the three solutions have self- acceleration. There are a finite density for GBIG1 
and GBIG2 but GBIG3 has big bang singularity. For negative non-minimal coupling, 
these solutions are not physically reliable since for all of them the density is negative. 



70 




100 



energy density 



Figure 7: Solutions of the Friedmann equation with a positive tension brane in a Minkowski 
bulk. 7 = 0.003, a = 10, t/j = 0, ol(4>) = 0.2 and —0.2 for dot and solid curves respectively. 

Note also that <7j )e and o\ are quantities for which /Zj je = and fii = respectively. 
From equation (37), <7j je are given in terms of 7 and a(0) by 

2a 3 (0)-9a(0)7 + 18a(0)^h±2(« 2 (0)-3 7 + 6 7 ^«:l) f 
= 27^ ' (40) 

and from equation (38), o\ = 0. In the non-minimal case, the maximum value of 7 is 
Imax = „ fdifl) 2 - For 7 = f d ifl 2 , GBIG1 branch disappears, since hi = h e at -f max . 
Actually the point hi = h e is now a point of inflection. The requirement to have real 
value for the square root in equation (36) and (37) leads us to the following relation 

7 < ~ „ rfaW 2 - (41) 

Based on this relation, the range of variation of 7 depends on the non-minimal coupling 
coefficient directly. The latest observational constraints on the values of 7 are listed in 



17 



Table 1: The values of 7 in the different test. 



Test 


7 


SNIa 
SNIa+LSS 
SNIa+LSS+H(z) 


0.0278lg : g?g 

o.ooo±8:8SB 
o.ooo±8:8S 



table 1 [55]. In this framework, we can constraint non- minimal coupling of scalar field 
and induced Ricci scalar using observational data. Considering a conformal coupling of 
the scalar field and induced gravity on the brane defined as a(<j>) = (1 — £0 2 ), we can 
constraint £ based on the constraints imposed on 7 presented in table 1. A detailed study 
of constraints on non-minimal coupling of scalar field and gravity based on various obser- 
vational data and theoretical techniques are summarized in reference [38] ( see also [56-64] 
for more details). For a time-independent scalar field, a((f>) will be time-independent also. 
Using the result of SNIa+LSS+H(z)( the third line of table 1) for 7 and relation (41), we 
obtain a > 0.01 and a < —0.01. Since a(<f) = (1 — £0 2 ), we find 




for a > 0.01 and 

> ]j and <fi < — 

for a < —0.01 respectively. The values of £ are constraint to be in the range of £ < O.990 -2 
and £ > 1.O10 -2 . As we see these constraints are dependent on the scalar field and this 
is reasonable since dynamics of scalar field essentially affects its coupling to gravity. In 
this viewpoint, variation of gravitational coupling as a field in scalar-tensor gravity can 
be attributed to variation of non-minimal coupling. 

4.3 AdS Bulk (ip ^ 0) with Brane Tension((j ^ 0) 

For ip ^ 0, the bulk is AdS since A 5 7^ 0. We should solve equation (23) for this case. 
The condition h — gives two solutions 

Hb,c = T^P(l + \)-^ (42) 
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where fi c is the density at the point that GBIG3 collapses (corresponding to plus sign), 
while fib is the density of the new bouncing point for GBIG4 (corresponding to minus 
sign). In the previous subsection (Minkowski bulk) there is a loitering point. Here this 
point separates into bouncing and collapsing points. It is interesting that when ip = —5, 
the GBIG4 branch disappears. Although the exact value of ip for disappearance of this 
branch is not important and depends on the choice of parameters, the fact that in principle 
this branch can be disappeared in a suitable parameter space is an important result. To 
obtain turning points of the branches, we calculate ^ as follows 



d{h 2 ) 



l + 7(^ 2 + f) 3 7 (/i 



f) + l 



2(a(0)/i 2 -(fi + fi' + a) + 2*$&h4) 



a{<t>ftf + 3a(0)^fK + + *)(«(*) + + 2( 



da(4>) ^2\2 
dr 



K i 



ifM + // + a) + 2^h4) 



2(a(0)/i 2 

d< - p+fi ) — 0, these points can be obtained by solving the equation 



(43) 



0. 



By substituting 

There are three roots (which are presented in Appendix A ) two of which are complex. 
Using equation (23), for AdS bulk satisfies the following equation 



: \A 2 - ~i 



(44) 



r 



+ 



1 + 2aa - i/ry(l + f^7) ~ 4 ( 



^2 



r 



.da L ,2\ 



V>(1 + + <7 2 ' 



0. 



- - - - (^) 

In comparison with the minimal case [35], the behavior of the branches are changed 
considerably. To see these differences, we obtain numerical solutions of the above equation 
for different values of ip. The results of this calculations are shown in figures 8, 9 and 10 
respectively. 

As these figures show, a(4>) as non- minimal coupling of the scalar field and induced 
gravity on the brane, controls the initial density and the age of the universe in this sense 
that these quantities are sensitive to the proposed value of the non-minimal coupling. 
For large values of a((f>), the universe age and its initial density are higher than the case 
with a small value of a((j>). Moreover, by increasing a((j>), one of the solutions, that is, 
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Figure 8: The solutions of the Friedmann equation with a positive brane tension in AdS bulk. 
We have chosen 7 = 0.003, a = 10, V = -5, a(4>) = 0.2, 0.3. 
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Figure 9: The solutions of the Friedmann equation with a negative brane tension in AdS bulk 
and in two different scales. We have set 7 = 0.003, a = —4, ip = — 1, a(4>) = 0.2, 0.3. 
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Figure 10: The solutions of the Friedmann equation with a negative brane tension in AdS 
bulk and in two different scales. 7 = 0.003, a = —4, ip = —1, a{<p) = 1 i.e. minimal case. 

GBIG4 disappears from the set of the solutions (note that these results are obtained for 
constant values of tp and a while ct(<p) is variable). The GBIG4 branch gives a bouncing 
cosmological solution. A bouncing universe goes from an era of accelerated collapse to 
an expanding phase without displaying any singularity. In the bouncing universe, the 
equation of state parameter of the matter content, uj, must transit from uj < — 1 to uj > — 1 
[65] . However, current observational data show that the equation of state parameter uj was 
larger than —1 in the past and is less than —1 today [66,67]. In our framework, we have 
seen that in a suitable domain of parameters space, by increasing a((f>) values the GBIG4 
branch containing a bouncing cosmology will disappear. Even for a constant a((p) there is 
no bouncing solution for any values of ip. For example with ip = —5, this branch disappear 
completely. These arguments show that inclusion of non-minimal coupling of scalar field 
and induced gravity on the brane can be used to fine-tune braneworld cosmological models 
in the favor of observational data. 

5 Time Evolution of the Branches 

In this section we discuss more general case of a time-dependent non-minimal coupling. All 
arguments presented in the preceding sections can be reconsidered in this time-dependent 
framework. We focus on Minkowski bulk (ip = 0) with brane tension(<r 7^ 0) for instance. 
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Starting with Friedmann equation (35), we try the following ansatz 



(j>(t) oc t~ 



(46) 



in order to investigate the late-time behavior of this scenario. As has been mentioned at 
the end of the subsection 4.2, the values of £ are constraint to be in the range of £ < O.990 2 
and £ > 1.O10 2 . By adapting these conditions and choosing a(<f>) = 1 — £0 2 , equation 
(35) can be rewritten as follows 



where proportionality constant in (46) has been set equal to unity. The constraints on 
£ are now time- dependent as £ > 0.99t~ 2u and £ < LOU -2 ". The results of numerical 
solution of this equation are shown in figures 11 and 12 for different values of £. Note that 
three graphs of figure 11 ( and also figure 12 with a different value of non-minimal coupling 
coefficient £) are corresponding to three branches of figure 6, but now with time variation 
of non-minimal coupling. These solutions are various possibilities of GBIG scenario with 
a time varying non-minimally coupled scalar field on the brane. For instance, figures 11a 
and 12a are corresponding to GBIG2 branch of the scenario. On the other hand, figures 
116 and 126 are corresponding to GBIG1 branch and finally, 11c and 12c are corresponding 
to GBIG3 branch. The main point to stress here is the fact that in the presence of explicit 
time evolution of scalar field, these branches show more or less the same late-time behavior 
as discussed in previous sections. 



(1 + 7 /i 2 ) 2 /i 2 = (1 - £r 2 > 2 -(// + // + a - 2/c^t(- 2 "- 1 )/ l ) 



(47) 
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Figure 12: Solutions of the Friedmann equation with a negative tension brane in a Minkowski 
bulk. We have assumed 7 = 0.003, a = -4, ip = 0, v = 0.9, «§ = !,£ = 2^ 2 . 
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Finally the issue of stability of the self-accelerated solutions should be stressed here. 
It has been shown that the self-accelerating branch of the DGP model contains a ghost 
at the linearized level [68]. The ghost carries negative energy density and it leads to 
the instability of the spacetime. The presence of the ghost can be related to the infinite 
volume of the extra-dimension in DGP setup. When there are ghosts instabilities in 
self-accelerating branch, it is natural to ask what are the results of solutions decay. One 
possible answer to this question is as follows: since the normal branch solutions are ghost- 
free, one can think that the self-accelerating solutions may decay into the normal branch 
solutions. In fact for a given brane tension, the Hubble parameter in the self-accelerating 
universe is larger than that of the normal branch solutions. Then it is possible to have 
nucleation of bubbles of the normal branch in the environment of the self-accelerating 
branch solution. This is similar to the false vacuum decay in de Sitter space. However, 
there are arguments against this kind of reasoning which suggest that the self-accelerating 
branch does not decay into the normal branch by forming normal branch bubbles [68]. It 
was also shown that the introduction of Gauss-Bonnet term in the bulk does not help to 
overcome this problem [69]. In fact, it is still unclear what is the end state of the ghost 
instability in self- accelerated branch of DGP inspired setups (for more details see [68]). 
On the other hand, non-minimal coupling of scalar field and induced gravity provides a 
new degree of freedom which requires special fine tuning and this my provide a suitable 
basis to treat ghost instability. As we have shown, non-minimal coupling of scalar field 
and induced gravity has the capability to remove bouncing solutions. It seems that this 
additional degree of freedom has also the capability to provide the background for a more 
reliable solution to ghost instability. This issue deserves as a new research program. 



6 Summary and Conclusions 

DGP model modifies the IR sector of general relativity. In the UV limit of a reliable 
theory, stringy effects should play important role. In this viewpoint, to discuss both 
UV and IR limit of the scenario simultaneously, the DGP model is not sufficient alone 
and we should incorporate stringy effects via inclusion of the Gauss-Bonnet terms. The 
presence of GB term removes the big bang singularity, and the universe starts with an 
initial finite density. Non-minimal coupling of scalar field and induced gravity on the 
brane which is motivated from several compelling reasons, controls the value of the initial 
density in a finite big bang cosmology on the brane. This is not the only importance of 
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non-minimal coupling of scalar field and induced gravity; non-minimal coupling provides a 
mechanism for generating spontaneous symmetry breaking at Planck scale on the brane. 
In this respect, non-minimal coupling of scalar field and induced gravity on the brane 
itself is a high energy correction of the theory and it is natural to expect that this effect 
couples with stringy effects in Planck scale. Investigation of the late-time behavior of DGP 
scenario with GB and non-minimal coupling effects provides a framework for constraining 
non-minimal coupling using recent observational data. In our model, the values of non- 
minimal coupling are constraint so that the values that £ can attain are constraint to 
be in the range of £ < O.990~ 2 and £ > 1.O10 -2 . In our setup, these constraints are 
dependent on the scalar field dynamics and this is reasonable since essentially dynamics 
of scalar field affects its coupling to gravity. One of the main outcome of our analysis 
is the implication of non-minimal coupling on bouncing cosmologies. In the bouncing 
universe, the equation of state parameter of the matter content, u, must transit from 
ou < — 1 to u; > — 1. However, current observational data show that the equation of state 
parameter u was larger than —1 in the past and is less than —1 today. In our framework, 
we have seen that with a suitable parameter space, by increasing a(<p) values, the GBIG4 
branch containing a bouncing cosmology will disappear. Even for a constant a((f>) there 
is no bouncing solution for any values of ip. These arguments show that inclusion of non- 
minimal coupling of scalar field and induced gravity on the brane can be used to fine-tune 
braneworld cosmological models in the favor of observational data. Although most of the 
arguments in the paper are based on a time independent non-minimal coupling, but as 
we have shown, inclusion of an explicit time dependence of non-minimal coupling will not 
change the physical nature of the solutions. Finally we have discussed the issue of ghost 
instabilities in self-accelerated solutions and possible impacts of Gauss-Bonnet term and 
non-minimal coupling on this issue. 
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APPENDIX A: Three Roots of Equation (44) 



hl = W _ ^ " ( 6(37 " «)( 2 - 3 ^ + 2 ^) 
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4(6(37-a)(2-37^+2a^)-4A 2 ) 3 + (-2167A-16207 2 ^+72aA+972 7 ^tt^-144^a 2 A 

' 2I67A - 16207 2 ^A + 72aA + 972^aA - U4ipa 2 A 



+ 



+ 16A 3 ) 2 }' , 

J / J 6 x 21(37 -a) 

+ 16 A 3 + [4 (6 (37 - a) (2 - 37^ + 2a^) - 4A 2 ) 3 + (-2 167 A - I62O7 V A + 72a A + 972 7 ^a A 



3n2 



-144^a 2 A + 16A 3 ) 



1 \ 3 
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3(3 7 



- ^(l +iv / 3)(6(37-a)(2-37^ + 2a^) -4A 2 )^J/ 
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